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Abstract 


We study the significance and characteristics of interfacial dynamical charge transfer at 
metal-organic interfaces for the organic semiconductor model system 1,4,5,8-naphthalene- 
tetracarboxylic dianhydride (NTCDA) on Ag(lll) quantitatively. We combine infrared ab¬ 
sorption spectroscopy and dispersion-corrected density functional theory calculations to ana¬ 
lyze dynamic dipole moments and electron-vibron coupling at the interface. We demonstrate 
that interfacial dynamical charge transfer is the dominant cause of infrared activity in these 
systems and that it correlates with results from partial charge and density of states analy¬ 
sis. Nuclear motion generates an additional dynamic dipole moment but represents a minor 
effect except for modes with significant out-of-plane amplitudes. 

Metal-organic interfaces are essential elements influencing the performance of molecular elec¬ 
tronics applications like light-emitting devices or field-effect transistors. Here, electronic exci¬ 
tations via charge-transfer excitons often play a key role. The formation mechanisms of these 
excitons via electron-hole pair (EHP) formation at the interface thus need to be understood. 
The coupling of EHP to adsorbate vibrations (electron-vibron coupling) is one way to induce an 
interfacial dynamical charge transfer (IDCT), which delivers key information about the electronic 
structure at the interface. IDCT induced by an orbital dipping in and out of the Fermi sea has 
been suggested for small molecules on metal surfaces.®^ This model has been applied to larger 
organic adsorbates on metal surfaces more recently.!^ The effect of such coupling on the vibra¬ 
tional line shape was analyzed by Langreth.^^ Despite its importance, there is no theoretical or 
experimental proof for IDCT as such and the relative importance compared to nuclear motion at 
a metal-organic interface beyond heuristic models up to now. We quantify this contribution by 
first-principles analysis and provide a rationale for the amount of IDCT for specific vibrational 
modes at the interface. 

Prerequisites for the occurrence of IDCT are (i) partial occupation of an adsorbate molecular 
orbital upon adsorption and (ii) strong electron-vibron coupling for relevant substrate electronic 
levels.l^ This was met in early studies for CO on Cu(lOO)^ and O2 on Pt(lll).^ Later, IDCT 
was observed in studies on fullerene IDCT, in conjunction with EHP excitation is also 

held responsible for the strong damping of high frequency vibrations of adsorbates via energy 
dissipation, which was one of the driving forces for its analysis. PF°F^ I Line-width analysis for IR 
measurements then leads to an estimate of vibrational lifetimes.^ The occurrence of electron- 
vibron coupling is essentially a breakdown of the Born-Oppenheimer approximation. For systems 
with strong non-adiabaticity, IDCT thus leads to asymmetric line shapes in vibrational spectra.I^^^m 
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While a theoretical treatment of this asymmetry requires a time-dependent analysis,^ which is 
not feasible for the system investigated here, an empirical evaluation for NTCDA on Ag(lll) has 
been presented recently by Braatz et al.® 

In this Letter, we will present a microscopic view on the vibrational modes of 1,4,5,8- 
naphthalenetetracarboxylic dianhydride (NTCDA) on Ag(lll) based on experiment and compu¬ 
tation focusing on the contribution of IDCT to the observed infrared (IR) absorption intensities. 
We correlate this to results of charge-transfer analysis associated with the electron-vibron cou¬ 
pling. We demonstrate that IDCT is the dominating effect for electron-vibron coupling in this 
system, quantify the charge-transfer associated and verify that IDCT is crucial for symmetric 
ag modes to become IR-active, in stark contrast to the free molecule. The system chosen here 
constitutes a model system for a more general class of metal-organic interfaces with noble metal 
substrates and planar vr-conjugated organic semiconductor species and we expect the results to 
be transferable to systems with similar energy level matching. Quantification of electron-vibron 
coupling has been identified as a major challenge for these interfaces.!^ 

Experimentally, infrared absorption spectroscopy (IRAS) is the method of choice to extract 
information regarding IDCT since it provides accurate vibrational energy and intensity measure¬ 
ments for modes originating at the interface. Furthermore, IRAS adheres to strict selection rules 
and probes molecular layers in a non-destructive way. The IR intensity of a mode is proportional 
to the square of the change of dipole moment of this vibration with respect to the equilibrium 
structure (dynamic dipole moment) 

^IR « Adyn) (1) 

where /Xdyn results from contributions of nuclear motion (ahuci) srid dynamic charge transfer 
across the adsorbate-substrate interface (aidct) 

l^dyn ^‘IDCT T f^nuch (2) 

as depicted in Scheme [T]:. Upon adsorption, charge q is being transferred from the surface to 
the adsorbate. Changes in q along vibrational coordinate Q (dynamic charge transfer Aq) can 
be determined by partial charge analysis; they are proportional to the IDCT contribution: 

A/r'cr oc Aq. (3) 

We will show that the contribution /Xnuci can be derived from equation]^ if /xidct is the major 
component. Additionally, the contribution from nuclear motion can be estimated directly from 

Aniici AAQ) (4) 
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where /xaq gives the dipole moment change of the adsorbate distorted along Q without the 
presence of the surface. For an adsorbate on a metal substrate only the z-component of fidyn 
contributes to the intensity due to dipole selection rules. Therefore, /inud is determined by out- 
of-plane displacements of the adsorbate atoms (Scheme [^). 



Scheme 1: Most stable adsorption geometry of NTCDA on a bridge position on Ag(lll) in top (a) and 
side (b) view; (c) contributions to dynamic dipole moment. 


Atomic and electronic ground state structures of metal-organic interfaces are well represented 
by dispersion-corrected density functional theory (DFT).®'^The description of vibrational spec¬ 
tra has also successfully been achieved in the past.^^lES 

NTCDA represents a model system with typical properties for an important class of organic 
semiconductors suitable for metal-organic interfaces: (i) Highly symmetric planar atomic struc¬ 
ture, (ii) delocalized vr-electron system, (iii) strong binding to metal surfaces with multiple bonding 
mechanisms. The totally symmetric ag modes are IR-inactive for the free molecule (/idyn = 0). 
According to IRAS the vibrational spectrum of NTCDA on Ag(lll) is strongly dominated by 
these in-plane vibrational modes, which suggests a notable contribution to /idyn from IDCT.^^ 
Near edge x-ray absorption fine structure measurements showed strong electron-vibron coupling 
in this system which renders it highly suitable for the quantitative determination of IDCT.i^ 

Our computational investigations applied dispersion-corrected DFT computations at the GGA 
level (PBE-D3)^2iH2Zlas outlined in the computational details section. The sub-monolayer structure 
investigated here differs from a genuine relaxed monolayer insofar as the NTCDA molecules are 
separated by one extra Ag atom row so that hydrogen bonding between neighbouring adsorbates 
can be neglected and all NTCDA occupy identical adsorption sites. NTCDA preferentially adsorbs 
with the central C=C bond occupying a bridge position on the Ag(lll) surface as shown in 
Scheme[I^. The detailed discussion of structural parameters, which agree well with experimental 
data of x-ray standing wave measurements,!^ is presented elsewhere.!^ The molecule is bonded 
by van-der-Waals interactions together with direct oxygen-silver chemical bonds (Scheme [^) 
leading to a strong adsorption and thus fulfilling the first prerequisite for IDCT.i^l The partial 
filling of the LUMO upon adsorption is a well-known phenomenon for NTCDA^^Sland thus the 
second requirement for IDCT is also satisfied. 
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In Figurethe experimental and computed IR spectra of NTCDA on Ag(lll) are displayed. 
The IRAS spectrum refers to 0.15 monolayers deposited at 28 K. The low temperature ensures 
that isolated NTCDA molecules are present and the formation of denser islands is avoided, in 
accordance with the arrangement used in our theoretical analysis. Comparison to previously 
derived spectra for the relaxed monolayer of NTCDA on Ag(lll)l^ yield moderate frequency 
shifts and only slight intensity variations, most probably due to some minor structural modification 
as a result of the different local environment. We obtain very good agreement of experiment 
and theory for vibrational energies and relative intensities. Thus, we can trust the computational 
spectrum to reproduce all the essential features and base the further examination on these data. 



Figure 1: Computed (top) and measured (bottom) IR spectrum for NTCDA on Ag(lll). Intensity scaling 
(x5, inset) was applied to signals below 1200 cm“^ in the computed spectrum. Flags indicate symmetry 
of respective modes applying free molecule nomenclature. 


The determination of mode symmetry for the vibrations of the adsorbate has been achieved 
by comparison of the displacement patterns to the free NTCDA molecule. Although the molecule 
bends upon adsorption (Scheme [^) and the point group symmetry reduces from D 2 h to approxi¬ 
mately C 2 v, the respective symmetry labels of modes are maintained as this way of labeling is the 
most helpful and intuitive one. In Figure[^it is apparent that most of the intense IR bands belong 
to the totally symmetric ag irreducible representation. Notably, these modes are IR-inactive in 
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the free molecule since they do not exhibit dynamic dipole moments. 




Figure 2: Equilibrium Structure (ES) of NTCDA with positive (ES+AQ, orange) and negative (ES-AQ, 
light blue) distortion along mode 13 (1565.6 cm”^, ag, top) and mode 12 (1509.8 cm“^, b 2 u, bottom) 
together with the corresponding pDOS stemming from pz states of NTCDA. The atomic displacements 
have been amplified. 


Next, we analyze the contribution of IDCT to Iir by displacing the adsorbate along normal 
mode coordinate Q. For the most intense mode at 1565.6 cm“^ (ag) the displacement pattern 
and the corresponding pDOS of the vr-orbitals (composed of atomic orbitals) of NTCDA are 
displayed in Figure]^ (top). Shown is the pDOS for the equilibrium structure (black line) together 
with positively (orange line) and negatively (blue line) displaced structures. The partial occupation 
of the bands stemming from vr-orbitals is unveiled (Fermi level crossing the pDOS lines), thus 
reflecting the known partial LUMO occupation upon adsorption. For the displaced structures, we 
see a strong shift of the pDOS maximum and a respective change in density at the Fermi level. 
Thereby, we observe electron-vibron coupling in this system and identify a strong dependence of 
the electronic structure on the vibrational motion - a direct prove for IDCT derived from first 
principles calculations rather than heuristic assumptions. Analogous behavior is observed for the 
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other Eg modes in the adsorbate, albeit at a reduced extent. Analyzing a larger entity of modes 
(Table [^, we find that IR-active modes exhibit notable values for f^dyn as expected from equation 
[I] The high values for Aq for most IR-active modes are a strong indication that IDCT is the 
leading effect in the system investigated. 

Further support for this hypothesis can be found by analyzing the modes exhibiting different 
symmetries. Modes 3 (717.9 cm“^, bsu), 5 (813.0 cm“^, bsu) and 15 (1628.9 cm“^, biu) exhibit 
significant, albeit much lower intensities in the experimental as well as the computed spectrum 
(Figure Table in line with their lower Aq values. For the bsu-symmetric modes 3 and 5, 
we find a significant out-of-plane bending component as expected from symmetry considerations. 
This leads to a major contribution to ^dyn that can be quantified by determining /iaq (Scheme 
[^, which shows that mode 3 has the highest contribution from nuclear motion.l^ 

Although we now established qualitative relations between jjidyn and its contributions from 
IDCT and nuclear motion, we aim at a direct relationship to allow the quantitative comparison 
of both factors. Figure reveals that the charge transfer contribution Aq indeed strongly 
correlates with jj-dyn- No correlation is found with /iAQ- Although correlation does not imply 
causality, combined with the pDOS analysis it becomes clear that IDCT determines Hdyn of ag 
type of modes and thereby Iir. ^idct can now be derived from Aq by the linear correlation 
given in Figure]^ under the assumption that it is the leading term (/i/dct ~ Adyn)'^3rid is given 
in Table [!} From equation]^ the signed value for jj^nuci can then readily be determined as the 
difference of fXdyn and hidct- The sign convention is thus that a negative value indicates ^.nuci 
pointing in the opposite direction compared to Hidct- The good correlation of jJnud with the 
computed /xaq in Figurevalidates the approach. This linear correlation can now be employed 
to convert /iaq values into jj^^ud (^9- @) given in Table We can thus derive directly 
from /iAQ or as difference from Hidct, while the latter approach gives information regarding the 
relative sign of the contributions. 

For most modes investigated, we find /x/dct to be the leading contribution to /x^yn (Table [^. 
More specifically, for three out of the four modes with highest intensity (9, 10, 13), IDCT effects 
are an order of magnitude larger. The third-most intense mode 14 (1625.7 cm“^) also shows 
high IDCT contributions but lower Ijr. This can be understood by the opposite sign of the strong 
finud contribution, which stems from a substantial vibrational amplitude for the downward-bent 
acyl-oxygen atoms at the periphery of the adsorbate. Moreover, analysis of the displacement 
patterns unveils a mixed nature for modes 14 and 15 with dominant ag (biu) contributions. We 
suspect that this mixing of character is provoked by their close spectral vicinity. This also explains 
the surprising intensity of mode 15 since biu symmetry should render it dipole forbidden on a 
metal surface (/Xdyn along NTCDA long axis). This example shows the limit of the assignments 
based on molecular mode symmetries which is clear-cut for all other vibrations and the need to 
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study the contributions to /idyn in detail. Modes with medium intensity (6, 7, 8) and ag symmetry 
exhibit higher relative contributions from nuclear motion but in all cases IDCT is the leading term 
here as well, while mode 1 shows equivalent contributions from both mechanisms. Modes 3 and 
5 (bsu symmetry), exhibit opposite signs for both terms: This leads to a near-zero intensity for 
mode 5 {fXnuci ~ I^idct) snd high intensity for mode 3 due to nuclear motion {unuci > I^idct)- 


Table 1: Computed properties of vibrational modes. See Scheme 1 for definition of terms. 


No. 

~a 

sym.^ 

int.'^ 

/^dyn 

Aidct'^ 

Anucl'^ 

Aq® 


1 

648.1 

Eg 

0.004 

0.18 

0.08 

+0.11 

0.01 

+0.14 

2 

665.4 

b2g 

0.000 

0.00 

- 

- 

0.00 

0.00 

3 

717.9 

bsu 

0.058 

-0.59 

0.52 

-1.08 

0.11 

-1.14 

4 

720.8 

au 

0.000 

0.00 

- 

- 

0.00 

+0.01 

5 

813.0 

bsu 

0.002 

0.11 

0.22 

-0.18 

0.04 

-0.17 

6 

987.0 

ag 

0.002 

0.12 

0.13 

+0.06 

0.02 

+0.09 

7 

1104.2 

ag 

0.021 

0.37 

0.22 

+0.02 

0.04 

+0.04 

8 

1256.9 

ag 

0.038 

0.53 

0.35 

+0.05 

0.07 

+0.08 

9 

1345.5 

ag 

0.581 

1.95 

1.88 

+0.07 

0.42 

+0.10 

10 

1404.8 

ag 

0.133 

0.94 

0.92 

+0.09 

0.20 

+0.12 

11 

1435.8 

biu 

0.000 

0.02 

- 

- 

0.00 

0.00 

12 

1509.8 

b2u 

0.000 

0.00 

- 

- 

0.00 

0.00 

13 

1565.6 

ag 

1.000 

2.51 

2.44 

+0.17 

0.55 

+0.21 

14 

1625.7 

ag 

0.264 

1.27 

1.84 

-0.62 

0.41 

-0.64 

15 

1628.9 

biu 

0.024 

0.37 

0.57 

-0.19 

0.12 

-0.18 


a Vibrational modes in cm b Mode symmetries, c IR intensities normalized to highest value. 
d Dipole moments in Debye, e Charges in e. 


For modes without significant intensities (i.e. 2, 4, 11, 12 in Table [^, none of the two 
mechanisms is active. Exemplary, for the b 2 u-symmetric mode at 1509.8 cm“^, the pDOS is 
represented in Figure (bottom). No dynamical change of the electron density at the Fermi 
level is found in line with negligible contributions from hidct- Thus, although several of these 
vibrations are IR-active in the free molecule they do not show any intensity on the Ag(lll) 
substrate due to screening of dipoles parallel to the metal substrate. 

The relative contributions of nuclear and electronic effects found here agree well with the 
experimental observation regarding Fano line shapes of IR absorption bands which may display 
different signs of the asymmetry parameter.^ The main factor here is whether the dynamic dipole 
moments due to nuclear displacement and IDCT point in the same or opposite directions. Thus 
we can elucidate the interdependence of both mechanisms via a detailed analysis of the intensity 
data and derive the relative signs from ab initio computations. 






Figure 3: (a) Correlation of charge transfer (Aq) and fx^yn to determine //idctI the open rectangles refer 
to first order correction of /idyn due to nonzero /inuci- (b) Correlation of /iAQ with ^nuci to determine 

/^nucl- 


In conclusion, we have investigated electron-vibron coupling effects associated with vibrational 
excitations at metal-organic interfaces for the model system NTCDA on Ag(lll). Using infrared 
absorption spectroscopy and employing density functional theory based bonding and vibrational 
analysis we derived unequivocal evidence for the dominating role of interfacial dynamical charge 
transfer (IDCT) for dynamic dipole moments and associated infrared activities. Nuclear motion 
(out-of-plane bending) is found to be a secondary mechanism only. The magnitude of IDCT for 
totally symmetric vibrational modes can be understood by pDOS shifts and derived from partial 
charge analysis. Contributions from nuclear motion can be determined reliably if IDCT is the 
leading term. This provides a simple, quantitative measure for the complex phenomenon. Our 
general approach will facilitate the analysis of other relevant metal-organic interfaces with similar 
matching of energy levels in the future. On the basis of the thus obtained improved understanding, 
fine-tuning of the atomic and electronic structure for these types of interfaces becomes feasible. 


Methods 

Computational details 

DFT computations at the GGA level (PBE-D3)*^^^ were carried out in a plane wave approach 
using the projector augmented wave method with an energy cutoff of 350 eV and a T-centered (3 
3 1) Monkhorst-Pack k-space grid within the code VASP 5.2.12.^23H36|-p|^e adsorption structure of 
one NTCDA molecule on a four-layer Ag(lll) slab (lattice parameter determined as a = 4.073 A) 
was optimized in a 4 x 4 supercell. IR spectra for relaxed structures (without employing empirical 
scaling factors) were derived via finite-differences calculation of a partial Flessian matrix, displacing 
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adsorbate atoms (a negligible influence of substrate atom displacement was found) and computing 
intensities based on the z-component of /idyn- This level of approximation was successfully used 
before for this system.^ IDCT has been deduced by displacing the equilibrium structure along 
the respective normal coordinates in positive (+AQ) and negative (—AQ) direction. For the 
displaced structures, projected density of states (pDOS) were derived as well as atomic partial 
charges with the Natural Population Analysis scheme.l^^*^ Tests with the electron density-based 
Atoms In Molecules (AIM)I3S! scheme gave comparable results for Aq. 

Experimental details 

The IRAS experiments were carried out in a UHV chamber at a base pressure of 5 x 10“^^ mbar. 
The Fourier-transform infrared spectrometer (FTIR) was a Bruker IFS 66v/S with evacuable 
optics. The spectra were recorded using a liquid N 2 - cooled Mercury-Cadmium-Telluride detector 
(spectral range 600 — 5000cm“^). In order to reduce the noise level, the displayed spectrum 
represents the sum of several consecutive measurements obtained at a sample temperature of 
28 K and using polarized light; in total, 12500 scans at a spectral resolution of 2cm“^ have 
been co-added for the sample and the reference spectra, each. The Ag(lll) crystal was cleaned 
by Ar+ - sputtering (30 min at Uion = 700eV, lion = 3/iA, and Tgampie = 373K), followed by 
annealing to 773 K for 5 min. NTCDA was evaporated from a custom-made thermal evaporator 
at Tntcda = 390K (deposition rate of 0.2 ML/min) controlled by a PtlOOO temperature sensor. 
During evaporation, the background pressure typically increases by Ap = 1 x 10“^° mbar. 
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